A1 = readmatrix('NiI2_MT_COOL.rw.dat')


hold on
xlim([0, pi])
plot(x1,f1,'Linewidth',2,'Color',[0 0 0.8])
scatter(x1,y1,10,[0 0 0.8],'filled')
xlabel('Polarization')
ylabel('Photocurrent [pA]')
xlim([0,pi])
ylim([0.75,1.05])
xticks(0:pi/4:pi)
yticks(0.8:0.1:1)
xticklabels({})
yticklabels({})
%set(gca,'XMinorTick','on','YMinorTick','on')
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
hold off

set(gcf,'Position',[0, 0, 200, 200])